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We study the scattering of a soliton-like propagating particle with a wah of bound particles, in 
several strongly interacting one-dimensional lattice models with discrete degrees of freedom. We 
consider spin-polarized fermions (anisotropic Heisenberg spin chain), the fermionic Hubbard model, 
and the Bose Hubbard model, using precise numerical time dependent Density Matrix Renormal- 
ization Group techniques. We show that in all integrable models studied, there is no reflection. 
Instead, an incoming particle experiences particle- hole transmutation upon entry and exit of the 
wall, and travels inside the wall as a hole, analoguous to Klein tunneling, even though the dispersion 
is highly nonlinear and there is no external potential. Two particles are added to the wall on the 
incoming side and removed on the opposite side. For spin-polarized fermions a single transmitted 
particle thus shifts the wall by two lattice sites, in complete contrast to classical physics. For both 
Hubbard models, the wall shifts by one doubly occupied single site. In the nonintegrable models 
studied, the same process occurs in linear superposition with backscattering events. We demonstrate 
a corresponding fermionic quantum Newton's cradle and a metamaterial with "tachyonic" modes 
travelling faster than in an empty system. We present a possible atomic scale signal counter for 
spintronics. Our scenario should be realizable in future cold atom experiments. 



The investigation of time evolution in non-equilibrium 
situations is a fast-expanding frontier in quantum many- 
particle physics. With the development of relevant exper- 
imental techniques, e.g. in cold-atom setups [lj|2], situa- 
tions that were quite academic a decade or two ago have 
become accessible. The study of coherent non-dissipative 
dynamics far from equilibrium has received a particular 
boost. 

An emerging theme is the dynamics and interaction 
of excitations over simple states. For condensed-matter 
models known and intensively studied for many decades, 
surprisingly little is known about the dynamical evolu- 
tion of excitations that are simple but are not eigen- 
states. Obvious examples are the dynamics of a few in- 
teracting particles on lattices, or a few flipped spins in an 
otherwise ferromagnetic background. In Hubbard mod- 
els, wide interactions "repulsively bind" on-site pairs and 
other clusters or combinations of bosons [3 . In models 
with nearest neighbor couplings like spin chains, bound 
states ( "multi-magnons" ) exist even at small couplings 
and their non-equilibrium dynamics has been the subject 
of recent research [4]-[7] . Interacting models and dynamics 
with single site resolution are now within experimental 
reach [9]. 

Collision processes are fundamental for analysing the 
physical laws governing the dynamics of classical and 
quantum systems. A famous classical realization is New- 
ton's cradle, exhibiting the effects of energy and momen- 
tum conservation. A quantum version has been realized 
in cold atom experiments [2 . 

In this paper we analyse scattering in one-dimensional 
strongly interacting quantum mechanical lattice mod- 
els with discrete degrees of freedom using standard pre- 



cise techniques for coherent many body time evolution 
[TOl HI] . We find intriguing phenomena caused by the dis- 
crete quantum nature of the system. We consider three 
prominent models, spin-polarized fermions, the fermionic 
Hubbard model, and the Bose Hubbard model. We will 
first briefly describe the main phenomena, exhibited in 
Fig. [l] Then we analyse the physics behind the observed 
behavior, arguing that it follows from conservation laws 
and the discreteness of the models. Finally we discuss a 
fermionic quantum Newton's cradle, a metamaterial with 
a tachyonic mode, and other possible applications. The 
supplementary material contains details of the calcula- 
tions and further examples. 



Main phenomena. 

The setup of what we call quantum bowling consists of 
an almost stationary "wall" of particles sitting on consec- 
utive sites, on top of an empty lattice. The wall is hit by a 
single soliton-like particle (Fig.jl]). The ensuing dynamics 
are qualitatively independent of details of the initial state 
or magnitude of the couplings, within a large range. For 
integrable models [121 [13] (Fig. [l^,b), no backscattering 
occurs. Instead, the incoming particle is transmitted as 
a hole through the wall, i.e. there is particle hole trans- 
mutation. This situation resembles a manybody version 
of Klein tunneling [15 , even though the dispersion of 
contributing modes is strongly non-linear and there is no 
external potential. Strikingly, because of particle num- 
ber conservation, there are then two particles that remain 
stuck at the front of the wall, instead of the familiar clas- 
sical result as in Newton's cradle, where just the single 
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FIG. 1: Quantum bowling: an incoming single particle in 
Gaussian shape crosses a wall of bound particles on N con- 
secutive sites. It undergoes particle-hole transmutation and 
shifts the wall by 2 particles. The figures show particle density 
{ni(t)). (a) Spin-polarized fermions with Vjt — 10. The 10- 
site wall is moved by 2 lattice sites. Inset: A wall of just 2 sites 
is still moved (V/t = 20). (b) Fermi Hubbard model. Left: 
up-spin density, right: down-spin density. Wall of A/" = 10 
doubly occupied sites, IJ jt — 100. (c) Bose Hubbard model. 
Wall oi N —\% doubly occupied sites, U /t — 30. The incom- 
ing particle is magnon-like 11 . Lower part in (c): density on 
condition that a particle is present to the right of the wall. 



incoming particle would be added to the wall. When 
the transmitted hole exits the wall, it is converted into a 
particle again, so that two particles vanish from the right 
side of the wall. For spin-polarized fermions (Fig. [l^) a 
single incoming particle shifts the whole wall by two sites 
to the left, in contrast to the classically expected shift 
by a single site. The exiting particle itself is shifted to 
the right. At large coupling, the shift is by two sites and 
the shape of the wave packet is unchanged by the trans- 
mission. In both Hubbard models (Fig. [l]3,c), the wall 
consists of doubly occupied sites, and it is shifted by two 



particles, namely one doubly-occupied site. Again there 
is particle hole transmutation. In the Fermi-Hubbard 
model (Fig.jljD), an incoming up-spin fermion becomes a 
down-spin hole inside the wall. In non-integrable models 
like the Bose Hubbard model (Fig. [l]^), there is a finite 
probability for backscattering. However, the final state 
is essentially a linear superposition of (i) a backscatter- 
ing event and (ii) transmission with particle-hole trans- 
mutation and wall-shift as described above. Indeed, at 
large couplings there is very little entanglement between 
a transmitted particle and the shifted wall, even when 
the transmitted particle itself is in an entangled Gaus- 
sian state. In the following, we will discuss the models 
individually. 



Spin-polarized fermions. 

The dynamics of spin-polarized fermions is governed 
by the Hamiltonian 



tv 



t 



I i 

(1) 

where t and V/t parametrize the kinetic and interaction 
energy of the fermions on the lattice, respectively, c] and 
Ci are creation and annihilation operators at site i, and 
Hi = c\ci. This model is equivalent by a Jordan- Wigner 
transformation to the ID spin ^ Heisenberg XXZ model 
with anisotropy A = J^/J^; = V/2t. We use t = 1 as the 
unit of energy and inverse time for all models, and discuss 

V > 2t. We note that when starting from a product 
state, results are identical for either sign of U [16], and, 
by adapting the arguments in ref. [16], for either sign of 
V/t^ i.e. for attractive and repulsive models. 

The model is exactly solvable by Bethe ansatz. It con- 
tains eigenstates made up of strings of M bound particles 
in an otherwise empty lattice [IKSJI^, which become more 
compact as V/t increases. Their energy is [12] 

E = 2t(cosh(M0) - cos K)smh ^/sinh(M(/)) (2) 

where coshcj) = V/2t and K is crystal momentum. Their 
time evolution is dominated by [5 a maximum veloc- 
ity 2t sinh0/sinh(M^) which decreases like ^yj^^M-i for 

V ^ 2t. At large V/t and large M, this velocity becomes 
exponentially small; M-strings are therefore stable on a 
very long time scale even when prepared at a fixed posi- 
tion in space. 

In the initial state we prepare a "wall" , a product state 
of A/" = 10 consecutive particles on an empty lattice. This 
is not an eigenstate, but it is instead made up of all sub- 
divisions of A^ into M-strings. Since smaller strings have 
larger velocities, they will in time "evaporate" [8^ from 
the wall. At large V/t the dominant contribution [7] is 
an almost stationary M = N string, followed by the sub- 
division into M = 9 and a single fermion with maximum 
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FIG. 2: Small anisotropy: density of spin-polarized fermions 
at y/t = 2.2 (A = 1.1 in the equivalent Heisenberg model); 
evaporatively cooled initial state. Inset: scattering region at 
larger scale. 



velocity 2t [TT] . In the initial state we also prepare a sin- 
gle particle in a Gaussian superposition pJJ of width 4 
and momentum — 7r/2 (or alternatively as a local particle 
at the left boundary [11 ) which moves towards the wall 
with velocity 2. 

Being integrable, the tV model contains a macroscopic 
number of conserved quantities, of which particle num- 
ber n = ^lUi^ energy and energy current (thermal 
current) = ^ijf are especially relevant, with 



3i 



itV{cl_^ci{ni^i - -) + {m 



+ it'^cl_^ci^i 



1. 
h.c. 



)c|q+i) 



Since the interaction is local, the energy current is con- 
served separately in spatially disconnected regions. 

The initial state has a finite right-moving energy cur- 
rent, localized around the incoming particle [11 . This 
current is conserved after the particle contacts the wall. 
Yet the wall is "full" and cannot accomodate an addi- 
tional particle. Only a hole can move through the wall. 
Inside the wall the hole possesses the same energy cur- 
rent and the same energy as the incoming particle on the 
empty lattice because of particle-hole symmetry. How- 
ever, because of particle number conservation, a particle- 
hole pair must have been created, and there are now two 
particles located to the left of the wall. Because of cur- 
rent conservation and energy conservation, they cannot 
move: if there was a backscattered particle, it would have 
additional energy and a left moving energy current, which 
would have to be compensated by a right moving current 
from another right moving hole, which would have fur- 
ther additional energy, so that energy conservation would 
be violated. 

When the hole exits the wall, it has to become a parti- 
cle again, so that two particles are taken away from the 
right side of the wall. Overall, the wall moves by two sites 
to the left in this quantum mechanical process, contrary 
to the classical situation. 

The particle-hole creation resembles Klein tunneling 
[15]. However, here we have no external potential, but 
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FIG. 3: Density of spin-polarized fermions with next nearest 
neighbor interaction. Main figure: Nonintegrable model with 
V = 10, V2 = 1-2 and a 10 site wall. Bottom part: density 
on condition that a particle is present to the right of the wall. 
Inset: Integrable model eq|3]with 1/ = 10, a = 0.1 and same 
wall width. 



instead a many-body effect. In contrast to Klein tunnel- 
ing, the dispersion is not linear, but a cosine. When one 
starts with an initially localized incoming particle [TT], 
all momenta contribute, yet particle- hole transmutation 
and the overall features of the transmission process are 
the same. 



Small couplings; narrow walls 

Particle-hole transmutation and the shift by two sites 
in the tV model are dictated by conservation laws and the 
Pauli principle, not by especially strong coupling. Indeed, 
a coupling of e.g. V/t = 2.2 (A = 1.1) still exhibits the 
same effects. At such smaller V/t, M-string eigenstates 
are spatially more extended. An initially prepared wall of 
N consecutive particles evaporates more than at larger 
V/t and becomes wider. Yet Fig. [2] shows that lack of 
backscattering and shift of the wall by 2 sites still occur 
at V/t = 2.2. 

Amazingly, a wall of only N = 2 sites already shows 
the same phenomena, including a shift by two sites, i.e. 
by the full wall thickness (Fig. [l^, inset). 

Role of integrability 

Conservation laws are essential for the observed effect. 
We probe the role of integrability by studying a nonin- 
tegrable model = Htv + "^2 (Fig. [3|, in 
which the energy current is not conserved. Now there is 
indeed backscattering. One might suspect the presence of 
a next-nearest-neighbor coupling to be responsible. How- 
ever, when one takes Htv and adds the conserved thermal 
current with next-nearest-neighbor terms. 



— Hi 



tV ^OLJ 



(3) 



another integrable model results, which does not show 
backscattering (Fig. [sj inset). This provides strong ev- 
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FIG. 4: Bipartite entanglement entropy Sab for the scatter- 
ing in Fig. [T^. Right inset: semiclassical picture of the time 
evolution of particles, demonstrating both left-shift of the wall 
and right-shift of the travelling particle. Space and time di- 
rections are as in the main figure. Left inset: spatial density 
distribution at large x; without wall (shifted right by 2 sites), 
and after passage through a wall. 



idence that integrability is indeed closely connected to 
the observed lack of backscattering. 

It remains an open question whether conservation of 
n^E^j^ and the restricted local Hilbert space are suffi- 
cient to suppress backscattering, or whether full integra- 
bility is necessary. It would be very interesting to study 
a nonintegrable model which conserves the above quan- 
tities, if such a model exists[T7j. 



Quantum mechanical nature of final state 

In the bottom part of Fig. [TJc) and Fig. [sj the wave 
function has been projected (and then normalized) onto 
Fock states in which exactly one particle is present to 
the right of the wall, i.e. onto the case that the incoming 
particle was transmitted. Then no reflection is visible in 
the reflected component and we see that now the wall is 
shifted by 1 (resp. 2) site. 

Further insight is gained from the bipartite entangle- 
ment entropy ^8|, Sab = — tr(pA logp^), where pA = 
tr Bp, and p = is the total density matrix. Sab 

quantifies the number and strength of linear superposi- 
tions between A and B. When it is zero, then is 
a product state \iP)a\'4^)b' In Fig. |4]we show Sab as a 
function of time and of the position of the subsystem cut. 
It is dominated by the entanglement inside the travelling 
Gaussian particle. Additionally, the slight evaporation 
of the wall is visible as light blue cones emanating from 
the wall boundaries. Strikingly, the amount of entangle- 
ment between the transmitted particle and the wall is 
hardly larger than on the left hand side of the wall [19]. 
Thus there is no, or very little, entanglement due to the 
outgoing particle. The outgoing particle, itself in an en- 
tangled Gaussian state, is therefore to good precision in 
a product state with the shifted wall. 

One can understand further details from a semiclassi- 
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FIG. 5: A fermionic Quantum Newton's Cradle: density of 
spin polarized fermions in a linear electric field at V — 10, 
performing Bloch oscillations. Initial Gaussian particle. V — 
10, 10-site wall. 



cal picture (Fig. |4j inset) in which the incoming particle 
is thought of as a single occupied site. Because of energy 
conservation, the closest that the particle can come to the 
wall is to a distance of one site; then a particle from inside 
the wall has to move to the left^ effectively propagating 
a hole inside the wall to the right. This picture implies 
that the propagating signal should experience a shift for- 
ward by 1 lattice site both upon entry and upon exit of 
the wall, thus overall the transmitted particle should be 
shifted by two sites in forward direction. Such a shift is 



4] and in the energy 
shows that at large 



indeed visible in the entropy in Fig. 
current [11 . The left inset in Fig. H 
V/t^ the Gaussian signal is moved forward by exactly 2 
sites without noticeable change of shape. At small V/t, 
we observe a smaller shift of about 1.3 sites, and we note 
that there will have been non-negligible additional scat- 
tering events with evaporated particles, which effectively 
widen the signal. 



Fermi Hubbard model 



The ID fermionic Hubbard model is specified by 
^fH = (4aC(i+i)a + h.c)j + ^ nitni; (4) 

1(7 i 

where U is the onsite interaction and a labels spin. It 
is widely used as a basic model of strongly correlated 
matter. Its coherent dynamics have recently started to 
become accessible in cold atom experiments [16 . 

The initial state in Fig. [ijb) contains = 10 consec- 
utive doubly occupied sites. Individually, they are re- 
pulsively or attractively bound by large \U\. They are 
however not mutually bound to each other. The outer- 
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most sites can therefore decay more easily than in model 
([T]), and |[/| > 30 is needed to clearly see the shift over 
background [11]. 

The model is integrable and Fig. [l] shows no backscat- 
tering at all. There is a conserved current [12 (slightly 
different from the energy current). Again, since the wall 
is already doubly occupied, only a hole can move through. 
Due to energy conservation, the number of doubly occu- 
pied sites has to be conserved. Therefore a hole of op- 
posite spin has to move. Since there is no direct nearest 
neighbor interaction, in a semiclassical picture the im- 
pinging particle can move up to the wall. Therefore, un- 
like the spinless fermion case, there is no forward jump. 
Indeed, we observe that transmission through the wall 
affects neither the shape nor the position of the wave 
packet [11]. 

Bose Hubbard model 

The dynamics of the Bose Hubbard model 

HbH = E (^l^^+i + '^■^■) +UY.n,{n,-l). (5) 

i i 

is now widely realized in experiments with ultra cold 
atoms in optical lattices [T]. Because of Bose statistics, 
there is no restriction on local occupation numbers. 

The initial state in Fig. [TJc) contains = 10 consec- 
utive sites with two bosons each, which are attractively 
or repulsively bound on each site. For the Bose-Hubbard 
model, the stability of doubly occupied walls grows with 
the wall width, as well as with |[/|, because the pairs can 
bind to each other when they are on neighboring sites 

m- 

The Bose-Hubbard model is non-integrable. The scat- 
tering result is similar to Fig. [Sj with partial backscat- 
tering. Because of energy conservation, any transmitted 
incoming signal has to go through the wall as a hole, 
leaving behind an additional particle and thus one addi- 
tional doubly occupied single site to the left of the wall. 
The projected part at the bottom of Fig. ^ shows that 
for these transmitted particles, the wall is indeed again 
shifted by 1 site, which remains visible for \U\ > 10. Note 
that inside the wall, the velocity of the hole is twice as 
high as outside, due to the double occupancy inside the 
cluster which renormalizes the hopping. 

Applications 

We discuss some immediate applications, making use 
of the wall shifts, the clean nature of scattering in the 
integrable tV model of spin-polarized fermions (or spin 
1/2 Heisenberg chain), and the large stability of walls. 

Fermionic quantum Newton^s cradle on a lattice. The 
continuous space "Quantum Newton's Cradle" [2 is one 




FIG. 6: Atomic scale signal counter and shift register: each 
passing particle shifts the wall by 2 sites. (Density of spin- 
polarized fermions, V = 10, 10-site wall.) 



of the most famous experiments with cold atoms. Here 
we construct a rather distinct lattice fermionic version 
by placing the system into an electric field with constant 
gradient, adding 0.06 jn^ to eq. (IT|). An incoming 
Gaussian particle then experiences Blocn oscillations [21 , 
whereas the wall is not affected noticeably due to its high 
mass. Fig. [5] shows the result: a periodic motion very 
similar to the classical Newton's cradle, except that at 
each impact, the wall moves by two sites instead of one. 

Qubits and atomic scale shift register. When several 
individual particles hit a wall in succession, the shifts 
add up quasi-classically, as shown in Fig. [6) In effect, 
the wall position counts the incoming particles, of po- 
tential practical interest, e.g. in spintronics applications. 
Furthermore, when a bound pair of particles impinges on 
the wall (not shown), it is transmitted inside the wall as 
a hole pair, shifting the wall by 4 sites after transmission. 

One way to encode a qubit with M-strings is to assign 
|0) to an M-string in a certain position and |1) to a similar 
M-string in a different position, ideally non-overlapping. 
We note that for large values of V/t and M, this qubit 
will decay only on timescales exponentiall small versus 
the inverse hopping 1/t. 

A sequence of walls and empty space, possibly of vari- 
ous thicknesses, can be interpreted as a sequence of bits. 
A suitable quantum-mechanical superposition of such se- 
quences can be seen as a sequence of qubits. An incom- 
ing single particle would shift the complete sequence of 
qubits coherently by 2 sites, without becoming entangled, 
making for a coherent qubit shift register [22] , 

Metamaterial, tachyon. We have shown above that 
the transmitted particle is shifted forward by 2 sites 
(Fig. |4j inset). When the particle moves through several 
walls in succession, as shown in Fig. [7| the individual 
shifts of the travelling particle add up. It moves with 
an average velocity higher than on the empty lattice. If 
one regards the latter as a "vacuum" with velocity vq, 
then the sequential walls act like a metamaterial with a 
tachyonic mode of velocity vod/ (d— 2), where d is the wall 
spacing {d = 7 in Fig. [t]). The results in Fig. [t] match 
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FIG. 7: Metamaterial with "tachyonic" mode: transition of 
a spin-polarized fermion through a set of 12 equally spaced 
4-particle clusters with an intra-cluster distance of 3 sites, at 
V/t—2^. Particle density shown. The forward shift of the 
travelling particle at each wall leads to an average velocity 
larger than the velocity of a particle on an empty lattice. 

this expectation precisely. 

Conclusions 

We have shown that the quantum mechanical trans- 
mission of a particle through a wall of neighboring 
particles exhibits surprising effects, namely pair creation 
with particle-hole transmutation and a shift of the wall. 
In the spin-polarized fermion or Heisenberg case, the wall 
shifts by two sites, and the transmitted particle jumps 
forward by two sites at large couplings. In addition, and 
independently, we find that there is no backscattering in 
the integrable models studied. These effects are due to 
conservation laws and the discrete nature of particles. 
They are therefore robust and, e.g., still occur at small 
anisotropics A > 1 in the Heisenberg model and for very 
narrow walls. The final state is close to a product state 
of a shifted wall and a transmitted particle of unchanged 
shape. Applications for spintronics may be possible, like 
an atomic scale signal counter or coherent shift of qubits. 
Last, but not least, the phenomena discussed should 
come within reach experimentally with cold atoms in 
optical lattices in the foreseeable future. 
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Supplementary material 

Methods 

We calculated the full quantum manybody time evo- 
lution of all models by time evolved block decimation 
(TEBD) [1 , which is closely related to the time depen- 
dent Density Matrix Renormalization Group (tDMRG) 
[2] . These techniques are based on Matrix Product State 
(MPS) representations and allow for high precision [3], 
especially with the relatively low entanglement in our 
simulations. We made use of local particle number con- 
servation and employed matrix dimensions of 200 - 400. 
We verified results against full diagonalization, where ap- 
plicable. For the Bose-Hubbard model, we limited the 
local occupation number on each site to Umax = 5, which 
was high enough to ensure that the results did not de- 
pend on Umax' 

Initial state. 

We started with an empty chain (length as in the fig- 
ures) with open boundary conditions. On this lattice we 
put a product state of TV consecutive singly or doubly 
filled sites to make the wall, by specifying MPS matrices 
accordingly. In addition, we started an incoming parti- 
cle, either in the shape of a Gaussian excitation, or of a 
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FIG. 81: Particle density for a single magnon-like excitation 
hitting a wall of A/" = 10 sites in the spin-polarized fermion 
model, V 10. 




FIG. S2: Energy current for a Gaussian particle hitting a wall 
of A/" = 10 sites (parameters of Fig.jTJi). 



localized magnon, as described below. 



Creation of Gaussian excitations. 
We followed ref. ^ and applied 

exp(--^^^-^^) exp{i{x - xo)ko) 4 

X 

^ J dkex.p{-— — ^ )exp(-zfexo)4 (S.l) 

to the state with the initial wall, with momentum ko = 
— 7r/2, width <j = 4, and center position xq as visible 
in the figures, and then normalized to create a single 
Gaussian shaped particle. Its velocity is |f Uo? where 
e{k) = — 2tcos(/c) is the single particle dispersion. With 
ko = — 7r/2, the resulting particle travels at velocity 2t 
with almost no dispersion, as seen in the figures. 



Single magnon excitation 

Alternatively, we started a "magnon" excitation by 
adding a single particle at the first left-hand site of the 
lattice. This procedure may be easier to implement 
experimentally. We also used it in the Bose-Hubbard 
model, because of better visibility of the transmitted hole 
inside the wall than with a Gaussian particle. 

The behavior of such an excitation is easiest to un- 
derstand by considering an empty lattice with a single 
particle in the middle at a site x = at time t = 
[5|. Since a single particle cannot interact, all models 
considered in this paper are equivalent to tight bind- 
ing fermions in this case. The initial state is \ip{0)) = 
|l).=o = Ek\^)k' Then m)) = eM-iHt)mO)) = 
exp{i2t^p cos php) \l)k = ^f^ex.p{i2tcosk)\l)k, 
which can be written as a Bessel function [5]. Each k 
mode travels with a velocity Vk = -2t^-^^ = 2tsmk < 
2t. Close to the maximum velocity v^ax = 2t, the most 
modes contribute, which produces a magnon, a linearly 
propagating wave distinctly visible in space-time figures. 
Modes which are further away from /c = 7r/2 produce 
additional oscillatory behavior. When starting in the 
middle of an empty lattice, a left-moving and a right- 
moving branch ensue, each with an overall probability 
of 1/2 of containing the particle. When started close to 
an open boundary, the particle travels only in one di- 
rection, away from the boundary, with unit probability. 
We note that the distinctly visible fastest branch has a 
velocity Vmax = 2t at its fastest edge, whereas the loca- 
tion of maximum intensity slowly moves away from the 
edge sublinearly. After finite times, the average velocity 
of the peak is therefore slightly smaller than Vmax^ but it 
converges towards Vmax for large times. 
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FIG. S3: Particle density for evaporative cooling of a 10-site 
wall at V — 2.2. The vertical scale has been cut off at n = 0.2 
for better visibility. 
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FIG. S6: Fermi Hubbard model. The spatial density distri- 
bution of the transmitted particle is independent of [/, and 
independent of whether there is a wall or not. 
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FIG. S4: Particle density at V = 10. Same as Fig. [TJi in 
the main paper, but with vertical scale cut off at very low 
n = 0.02 to make evaporation of wall visible. 



In Fig. SI we show the analogue of Fig. [T^ in the main 
text, with an initial magnon impinging on the wall, in- 
stead of a Gaussian particle. We see that the total in- 
tensity of the incoming particle now arrives over time in 
several waves. Each of them behaves similar to Fig. [l^. 
There is no backscattering, and the probability distribu- 
tion of the wall converges to a complete shift of 2 sites. 



Energy current for spin polarized fermions 

Fig. [S2] shows the energy current corresponding to 
Fig. [T^ in the main text. The current is conserved glob- 
ally, and in spatially disconnected regions also locally. 
The local nature of the incoming Gaussian particle and 
its jump forward at both edges of the wall are clearly 
visible. 



Evaporation and evaporative cooling. 

At V/t = 2.2 we prepared a wall of TV = 10 particles 
and let it evolve for time 100/t, during which particles 
evaporated and the wall became wider, closer to a linear 
combination of M-string eigenstates with mostly large 
M [6 , in a kind of evaporative cooling. Fig. S3 clearly 



shows the evaporation of single particles and of slower 
M = 2 bound strings. We then cut the remaining state 
to a width of 110 sites and started a Gaussian particle, 
as visible in Fig. [2] in the main paper. 

For comparison. Fig. [S4] shows the much smaller scale 
of evaporation at V/t = 10 for the time evolution in 
Fig. [Tk of the main paper. 
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FIG. S5: Densities in the Fermi Hubbard model at U — 30, 
similar to Fig. [TJ) in the main text. Left: up-spin, right: 
down-spin. The wall widens over time, but a shift of 1 lattice 
site (2 particles) is visible in the spatial density distribution. 



Fig. [S5] displays scattering for the Fermi Hubbard 
model, similar to Fig. ^ of the main paper, but at 
smaller coupling U/t = 30. The wall decays faster in this 
case. Its shift is still visible in the spatial density dis- 
tribution before and after the scattering. Fig. [S6] shows 
that, within the precision of our data, the shape of the 
transmitted particle is independent of U. It is not in- 
fluenced by additional scattering events with evaporated 
particles. This indicates that the particle emerges prac- 
tically unchanged, and without phase shift, from each 
scattering for a wide range of U values. 
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Interaction inversion symmetry 

In the appendix of ref.0it was shown for the fermionic 
Hubbard model that the time evolution is invariant under 
the transformation U — /7, for observables invariant 
under a combined transformation of time reversal and 
TT-boost Cq Cq+Q^ which iuclude particle density, and 
for initial states which only aquire a phase under the 
combined transformation. We note that the proof applies 
to density-density interactions in general, including the 
t^-model and the bosonic Hubbard model, and to initial 
product states. For an initial state including a Gaussian, 
it applies when —ko is also transformed. All our 

results therefore apply to both attractive and repulsive 
models, i.e. they are invariant under U —U, resp. 
V ^ -V. 
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